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ABSTRACT 

The stochastic forcing theory of Frankignoul and 
Hasselmann, 1977 is modified to include a mixed layer model. 
This enables the examination of the interaction between 
stochastic heat flux or wind stress components and the 
annual period surface heat flux. The presence of the 
stochastic heat flux component causes the average sea 
surface temperature to be 0.75 C higher then it would be 
with only the annual period component. It also delays the 
time of maximum surface temperature, and it causes the 
average mixed layer depth to be ten meters shallower. A 
stochastic wind stress component applied to an annual heat 
flux cycle produces a smaller sea surface temperature vari- 
ance, but it results in a more rapidly deepening and deeper 
mixed layer than achieved with an annual cycle and constant 
wind stress. The stochastic forcing model shows that the 
climatology of the seasonal thermocline is dependent on 
nonlinear interactions between the annual cycle and 
stochastic forcing. 
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I. INTRODUCTION 



The object of this paper is to show the effect of 
stochastic surface heat flux forcing on the climatology of a 
model seasonal thermocline. It was theoretically demon- 
strated by Hasselmann (1976) that if "weather" forcing can 
be considered to be white noise, then the ocean acts as an 
integrator. The result is a red-shifted response spectrum. 
Frankignoul and Hasselmann (1977) developed a simple model 
to test this theory. Here we shall expand on their result 
by considering seasonal variations and mixed layer physics. 

Climate researchers have long been aware of the 
pronounced variability associated with climatic records of 
sea-surface temperature (Namias, 1959; Wyrtki, 1965). This 
temporal variability spans a spectral range from a low 
frequency limit of the order of ice ages (one cycle per 
100,000 years) to a high frequency limit on the order of one 
cycle per month. ^ Higher frequencies than one cycle per 
month are taken to be stochastic in this study, whereas mean 
quantities and climatologies will have seasonal or longer 
frequencies. An understanding of the origin of this diverse 
climatic variability is a primary goal of climate research. 

Attempts have been made to link climatic changes to 
various external factors such as solar activity, perturba- 
tions in the Earth's orbit, or the increased atmospheric 
turbidity associated with volcanic eruptions. The primary 
difficulty with these theories is that the input-response 
relationships are not sufficiently pronounced to be immedi- 
ately obvious upon inspection of the appropriate time 
series . 



'The limits follow the definitions adopted in GARP 
Publication 16, 1975. 
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Recent interest, both commercial and tactical, focuses 
on the higher frequency sea surface temperature (SST) 
perturbations. The 1983 El Nino (seasonal frequency phenom- 
enon) proved disastrous to the commercial fishing industry 
along the west coast of the U.S. and South America. 
Tactically, the strategic defense of the U.S. is linked to 
the proper utilization of resources in the vicinity of SST 
anomalies and other thermal structures which alter acous- 
tical propagation through the oceans. 

Large scale SST-weather relationships have been discov- 
ered through the statistical analysis of extended time 
series (Bjerknes, 1966; Namias, 1969). Numerical experi- 
ments have demonstrated statistically significant modifica- 
tions of weather patterns in the tropics by critically 
located SST anomalies. These anomalies have been studied 
through the use of General Circulation Models (Shukla, 1975; 
Rowntree, 1976) and/or Statistical Dynamical Models. 

Much of the theory for ocean mixed layer modeling 
depends on the validity of two hypotheses. The first is 
that vertical mixing in the ocean mixed layer and entrain- 
ment at the thermocline occur as a result of atmospheric 
forcing: the heat fluxes across the air-sea surface and the 
surface wind stress. The second hypothesis is that the 
mechanical energy budget holds the key to understanding the 
ocean mixed layer dynamics (Garwood, 1979). 

The stochastic forcing concept is another approach to 
understanding the ocean SST response to atmospheric forcing. 
Lemke (1977) applied the stochastic forcing concept to a 
Budyko- Sellers type model to extend the concept to a global 
stochastic energy model applicable for time scales from 10 
to 10,000 years. Frankignoul and Reynolds (1983) examined 
the statistical characteristics of SST anomalies in the 
presence of a mean current. Their studies indicate the 
importance of surface heat flux and wind mixing fluctuations 
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on the production of SST anomalies from the available North 
Pacific Experiment (NORPAX) data. 

Hasselmann (1976) proposed the use of a stochastic 
forcing model to study the ocean response to atmospheric 
forcing. He assumed that, because the time scale of the 
ocean response is so much greater than the time scale of the 
atmospheric forcing, the atmospheric forcing can be consid- 
ered to be represented as a white noise time series. Via 
this representation, he showed that the slow climatic vari- 
ability of the ocean can be interpreted as an integral 
response to the random excitation provided by "weather" 
disturbances . 

In part two of the series on stochastic climate models 
(Frankignoul and Hasselmann, 1977), the stochastic forcing 
concept was applied to a very simple SST model utilizing a 
fixed depth mixed layer, "copper plate" ocean model. This 
simplest model of the ocean mixed layer assumes a completely 
mixed, horizontally homogeneous layer of depth h and temper- 
ature T overlying a seasonal thermocline with prescribed 
temperature structure. If the effects of salinity and hori- 
zontal advection are ignored, the layer parameters T and h 
are governed by two prognostic equations of the form (Kraus 
and Turner, 1967; Denman, 1973); 




( 1 . 1 ) 




( 1 . 2 ) 
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Equation (1.1) is a heat balance equation, where f| denotes 
the sum of: (1) the total flux of latent and sensible heat 
through the air-sea interface (dependent on the air tempera- 
ture Ta, humidity q and wind speed U via standard bulk aero- 
dynamic formulas); (2) the short and long wave radiation R 
through the sea surface; and (3) the heat exchange with the 
layer below. Equation (1.2) is an empirical relation for 
the entrainment of fluid from below the mixed layer and is a 
function of wind speed and other variables. 

Frankignoul and Hasselmann (1977) simplify the model 
further by retaining only equation (1.1), so that the mixed 
layer depth remains constant. Furthermore, only the contri- 
bution (1) to ft is included. Heat exchange with the atmos- 
phere is accomplished by assuming that the air-sea 
temperature difference is proportional to the north-south 
velocity V, via the relation (Ta-T)=KV, where K=constant. 
Therefore, a wind from the north will cool the ocean, and 
conversely a wind from the south will warm the ocean. Given 
the above parameterizations , the rate of change of SST is, 



The following values were assigned to the parameters of 
equation (1.3): 

K=0.25 CC/m/s) 

Ch = 10*’* 

B=3.0 

/^* = 1.25 X 10 ^ gm cm ^ 
gm cm*^ 

c; = 0.24 cal gm*'(*C)" 

CX=0.96 cal gm*'(‘'C)" 
h= 25 m 




(1.3) 
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A mixed layer depth of 25 m was selected as representative 
of lower mid-latitudes (30*N) on the basis of a comparison 
of the observed SST cycle with predictions by a "copper 
plate" model, (a value of h=100 m is used for calculations 
at station India 59*N, 19"W). With K=0.25, the root mean 
square (rms) air- sea temperature difference generated by the 
model was 1.25 C (Frankignoul and Hasselmann, 1977). 

The atmospheric model (Frankignoul and Hasselmann, 1977) 
is a barotropic model^ of quasi-geos trophic turbulence 
developed by Kruse (1975) for numerical investigations of 
nonlinear transfer within large scale atmospheric motions in 
an equilibrium state. It provides a statistically 
stationary wind field from which heat and momentum fluxes 
are inferred using bulk aerodynamic formulae. In this model 
there is no direct feedback from the ocean to the atmos- 
phere . 

Kruse's model represents a non-divergent , one layer 
atmosphere on a y9 plane bounded by two latitudes, and it is 
zonally cyclic. The dependent variable is the stream func- 
tion which obeys the equation. 



^Similar models are discussed by Lilly (1972) and Rhines 
(1975). 




where 




(1.5) 
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In equations (1.4) and (L.5), x and y are eastward and 
northward cartesian coordinates, Uo is a constant velocity 
of 5 m/s and J denotes the Jacobian operator (Frankignoul 
and Hasselmann, 1977). The boundary conditions are: 



and 

? pL, ps 1,2, ... (1.7) 



where, Lx=Ly=8800 km. Despite the high degree of simplifi- 
cation, the atmospheric model provided a statistically 
stationary wind field. 

Frankignoul and Hasselmann (1977) constructed time 
series of the stochastic atmospheric forcing from equation 
(1.3) at each gridpoint and calculated the SST changes via 
integration. To draw attention to the evolution of low 
frequencies in the SST fluctuations, the SST time series in 
Figure 1.1 have been low-pass filtered. Despite the 
simplicity of the coupled system, the white noise spectrum 
of atmospheric forcing results in a red- shifted SST 
response. Without a stabilizing negative feedback, the 
total SST variance is shown to grow indefinitely with time. 
With a properly chosen negative feedback, the response 
becomes asymptotically stationary. 

The purpose of the present research is to explore the 
role of seasonality in the ocean response to stochastic 
forcing. It is hypothesized that the seasonal variability 
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Figure 1.1 Simulated SST variations (without and with 
feedback effect) and total heat flux at x=0, y=Ly/4. 

The heat flux time series is subsampled at 5-day 
intervals and the SST time series are low-passed, 
usine a quadratic Lanczos filter (cut-off frequency 8 
X lO^HzJ (From Frankignoul and Hasselmann, 1977). 

in the surface heat flux (Qo)» as modified by a random 
forcing component (Q'), will cause seasonal variability in 
the mean mixed layer depth (h) and in the variance of mixed 
layer depth (h'"*). Variability in h will in turn modify the 
mean temperature (T) and the temperature variance (!'■*). 
This study is part of an investigation of increasing 
complexity : 

1) To a constant wind stress, impose a seasonal heating/ 
cooling (Qo) cycle plus a stochastic heat flux component 

(Q’); 
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2) Retain the seasonal heating/cooling (Q^) cycle minus 
the stochastic component, but add a stochastic component 
to the wind stress; and 

3) Consider the combined effect of an annual heat flux 
cycle with both a stochastic wind stress and heat flux 
component. This step represents the most realistic situ- 
ation, but care must be exercised to insure the proper 
correlation between the wind stress and surface heat flux. 

This paper examines the ocean response to random atmospheric 

forcing as described in investigations (1) and (2) above. 
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II. THE STOCHASTIC MODEL 



Frankignoul and Hasselmann considered a simple climatic 
system consisting of a slab oceanic mixed layer model which 
was driven by fluxes of heat and momentum across the air- sea 
interface. The governing equation (1.3) represents a linear 
approach to the study of climate variability due to 
stochastic forcing. The stochastic model now described is 
not only nonlinear in its approach, but it is designed to 
simulate mixed layer evolution on daily to annual time 
scales in response to atmospheric forcing. 

A. OCEAN MIXED LAYER MODEL INCLUDING SEASONAL THERMOCLINE 
Structure 

A potentially significant improvement over previous 
studies of stochastic forcing lies in the added complexity 
of the ocean mixed layer model. The latest version of the 
Garwood (1977) one-dimensional bulk model of the mixed layer 
of the ocean was chosen for several reasons. First, the 
model was readily available and required only minor modifi- 
cations for the task. Second, the model is thoroughly 
tested^ and reproduces mid- latitude boundary layer dynamics. 
Third, the model produces solutions of the turbulent kinetic 
energy components and for the terms of the turbulent kinetic 
energy budget. Future studies will more directly study the 
details of stochastic forcing on the turbulence components. 



^It is realized that there are no exhaustive tests of a 
model, but a very successful study of Ocean Station Papa 
simulations performed in 1981 supports the dynamic consis- 
tency of the model (Garwood and Adamec, 1981). 
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The model explicitly parameterizes a number of mixed- 
layer processes ignored by other models. The ratio of the 
mixed layer depth to the Obukhov length scale determines the 
fraction of wind- generated turbulent kinetic energy avail- 
able for mixing. The viscous dissipation is dependent on a 
local boundary layer Rossby number. The model uses separate 
equations for vertical and horizontal turbulent kinetic 
energy (see Garwood (1977) and Garwood and Yun (1979) for 
further details on the mixed layer model). 

The Garwood model is used here to simulate the evolution 
of the ocean surface layer on daily to annual time scales in 
response to atmospheric forcing. The horizontal gradients 
and horizontal advection are neglected in this one- 
dimensional model, and the version used in this project 
neglects vertical diffusion within the thermocline. In 
addition, the salinity is held constant over the water 
column. The equation of state is reduced to an approximate 
proportionality between temperature and buoyancy. 

In the absence of a vertical advection velocity, the two 
governing equations for the model become: 




( 2 . 1 ) 



and 




( 2 . 2 ) 
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where : 

T-mixed layer temperature 

Q^-mean annual surface heat flux (positive downward) 
Q' -random component of the annual surface heat flux 
at- temperature jump across the thermocline 
/OW-density of water 

C^-specific heat at constant pressure of water 
h-mixed layer depth 
e - ent r ainment 

Q, is prescribed according to: 



Where the amplitude A is equal to 150 Watts/m*. Q' is 
generated by means of a uniform random number generator and 
has a maximum potential amplitude equal to A. Therefore, 
the possible daily fluctuation in Q’ is from plus to minus 
A. 

The mixed layer temperature in this model changes as a 
result of the applied surface heat flux and entrainment 
mixing of colder water from below. Entrainment mixing of 
colder water from below the layer raises the potential 
energy of the water column. This energy for mixing is 
supplied by the mixed layer turbulence. Bulk layer models 
such as this one determine the rate of entrainment (e) from 
the turbulent kinetic energy budget. In the Garwood model, 
the time scale for viscous dissipation of turbulent kinetic 
energy depends on the intensity of the turbulence, the layer 
depth and on the planetary rotation. The rate of entrain- 
ment depends on the vertical turbulent kinetic energy avail- 
able and the temperature jump aT at the base of the mixed 





(2.3) 
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layer. Entrainment becomes zero when the vertical turbulent 
kinetic energy approaches zero. Thus the equation for the 
rate of entrainment (e) becomes: 




( 2 . 4 ) 



here 

w'=vertical turbulent velocity. 

<E>=mean mixed layer turbulent kinetic energy. 

«=thermal expansion coefficient. 
g=gravitational acceleration, 
m =dimensionless constant. 

The mean mixed layer turbulent kinetic energy <^> is a func- 
tion of both the surface wind stress and the surface heat 
flux. The relative distribution of the turbulent components 
is particularly sensitive to the sign and magnitude of the 
surface buoyancy flux. Positive buoyancy flux produces 
<w’w’> directly and increases vertical turbulent kinetic 
energy and entrainment. Negative buoyancy flux inhibits 
vertical mixing. Entrainment and mixing of water from below 
occurs only when there is sufficient turbulent kinetic 
energy available in the mixed layer. In this case there is 
active entrainment and mixed layer evolution is given by the 
prognostic equation (2.2). 

Under conditions of insufficient turbulent kinetic 
energy for mixing, the mixed layer retreats immediately to a 
depth hr. This retreat depth is determined by assuming 
steady state conditions for the turbulent kinetic energy 
budget and letting <w’w’> approach zero. The retreat depth 
hr is given by: 
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(2.5) 




> 



where 

U* is the friction velocity at the water surface, 
m, and nij are dimensionless model constants. 

With the prescribed surface heat flux (Q=Qo‘*'Q' ) and the 



and the mixed layer temperature and depth from (2.1) to 
(2.5). 

1. Reduction to " Copper Plate " Form 

The Garwood (1977) model is easily converted to a 
simple slab or "copper plate" model as used by Frankignoul 
and Hasselmann (1977). This is accomplished by setting the 
mixed layer depth (h) equal to some constant value. The 

result is that equation (2.2) becomes: 



Specifically, this implies infinite thermal conductivity 
(strong mixing) between Z=-h and Z=0, but zero conductivity 
(no mixing) below this. Given this condition, equation 
(2.1) reduces to: 



Equation (2.7) is effectively the form of equation (1.3) 
that Frankignoul and Hasselmann used. 



friction velocity U*, the model calculates <E> and <w' w' > 




( 2 . 6 ) 




(2.7) 
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The primary difference between (2.7) and (1.3) lies 
in the determination of the heat transfer between the ocean 
and atmosphere. Frankignoul and Hasselmann (1977) parame- 
terize the heat transfer as a function of the wind direc- 
tion, while Garwood (1977) specifies the surface heat flux 
directly. However, both (2.7) and (1.3) will produce equiv- 
alent stochastic forcing. The advantage of (2.7) is that 
the rate of change of SST is now dependent upon both the 
annual cycle and the random forcing, and facilitates later 
application with a predictive equation for mixed layer 
depth . 

Reducing Garwood’s model to a simple slab model 
permitted the recreation of the results of Frankignoul and 
Hasselmann (1977) for the case without negative feedback. 
These results will be discussed further in Chapter III. 

B. ATMOSPHERIC FORCING INCLUDING AN ANNUAL CYCLE 

The ocean model is driven by an annual period sinusoidal 
surface heat flux (Q^), a stochastically varying heat flux 
(Q’ ) and a constant wind stress. For computational conven- 
ience a model year is taken to be 360 days. The model is 
initialized at ’’spring transition” corresponding approxi- 
mately to the vernal equinox with a mixed layer depth of 100 
meters, a mixed layer temperature of 6 C and an underlying 
homogeneous water column of 5 C. The magnitude of the 
annual period surface heat flux was chosen so that the 
surface temperature would rise approximately 8 C and the 
seasonal thermocline would shallow to approximately 50 
meters at the times of maximum SST and maximum downward 
surface heat flux respectively. The necessary amplitude of 
the annual heat flux to achieve this was 150 Watts/m^in the 
case of a 1 dyne/cm* wind stress. The salinity is specified 
as a constant 33 7oo and it has no dynamical significance in 
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this study. The initialization described approximates 
conditions at Ocean Station Papa at 50®N, 145‘’W. To allow 
the model to reach a statistical equilibrium, the model is 
actually integrated to 720 days, and the second model year 
is used for calculations and graphic display. 

The stochastic forcing is supplied by a random number 
generator. A known but different seed value is specified 
for each case, which allows duplication of the random 
forcing series. The stochastic forcing is applied to the 
surface heat flux such that the random forcing can on any 
day equal (positively or negatively) the amplitude of the 
annual surface heat flux cycle. 

C. CALCULATION OF STATISTICAL VARIABLES 

At this point it becomes necessary to define the nomenc- 
lature to be used in discussions of the model results. 
Model outputs are sea-surface temperature, mixed layer 
depth, the average temperature variance T'^ and average 
mixed layer depth variance h'^. The variables To and ho 
represent the ocean response in temperature and mixed layer 
depth to the annual period sinusoidal heat flux cycle 
without the stochastic component. The differences T-T® and 
h-ho, between the mean surface temperature or mixed layer 
depth with random forcing and the mean surface temperature 
or mixed layer depth without random forcing provide insight 
into the effects of the stochastic forcing theory. 

An overbar denotes an ensemble mean of the daily values 
of the variable divided by the number of cases summed. 
Thus, the output statistics represent daily averages over 
some climatological time frame. The variance of temperature 
and mixed layer depth can be computed from the following 
equations obtained after performing a Reynolds 
decomposition. 
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( 2 . 8 ) 






and 



= -(h)‘ 



( 2 . 9 ) 



Thus, the variances are the result of the difference between 
mean quantities. 
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III. THE STOCHASTIC MODEL EXPERIMENTS 



A. PROCEDURES 

1 . General Method 

The stochastic forcing is applied to 101 cases, 
which corresponds to 101 model years. The first case or 
model year represents the nonstochastic case that determines 
the daily value of T» and h<> • The next 100 cases contain a 
daily random forcing component to the annual surface heat 
flux cycle. From these cases, the averages and variances of 
the temperatures and mixed layer depths are calculated. The 
average temperature and mixed layer depth are computed daily 
via a simple summation over the 100 stochastically forced 
cases. For the variance computation see equations (2.8) and 
(2.9). 

2. Reproducing Frankignoul and Hasselmann* s Result 

Frankignoul and Hasselmann (1977) applied the 
stochastic forcing concept to a simple slab or "copper 
plate" ocean model as given by equation (1.3). The results 
obtained from the Garwood model as in equation (2.7) confirm 
the hypothesis of Hasselmann (1976) that the ocean acts as 
an integrator of the white noise atmospheric forcing and 
produces a red-shifted response spectrum, that is the 
spectral response of the model is shifted toward the lower 
frequencies, and appears geophysical in nature (Figure 1.1). 
Frankignoul and Hasselmann further showed that, without a 
stabilizing negative feedback, the total SST variance grows 
indefinitely with time, as suggested by Hasselmann (1976). 

Figure 3.1 displays the daily model response to the 
sinusoidal annual heat flux cycle without stochastic 
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Figure 3.1 Response to annual cycle, .a) 
net surface heat flux (positive upwards). b 
ocean response to the annual period surface 
without stochastic forcing; h„ equals lUU m. 
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forcing. It can be seen that the annual surface heat flux 
cycle begins at vernal equinox, and adds heat to the mixed 
layer with the maximum heating rate occurring at 90 days. 
Heat continues to be added to the mixed layer for an addi- 
tional 90 days, and at day 180 the cycle reverses and heat 
is withdrawn from the ocean, with the maximum withdrawal 
rate on day 270. Heat continues to be removed from the 
ocean until the cycle is complete on day 360. 

Although the one-dimensional model is designed so 
that the mixed layer temperature will rise approximately 8 C 
with realistic mixed layer depths, it is evident that fixing 
the mixed layer depth reduces the ocean temperature response 
to only approximately 3.5 C. Since the mixed layer depth is 
constant there is no interaction between the minimum depth 
at the time of maximum surface heat flux. Instead, the heat 
input is spread over the entire 100 m, and the temperature 
increase is only about half of the expected value. 

The effects of the stochastic forcing for a typical 
case can be seen in Figure 3.2. The SST would statistically 
be expected to be larger than T© one half of the time and 
smaller than T^ the other half. In this case the SST is 
slightly smaller than T„. 

Figure 3.3 shows the significant temperature statis- 
tics calculated from the model. Hasselmann (1976) and 
Frankignoul and Hasselmann (1977) showed that the SST vari- 
ance in the absence of negative feedback will grow indefi- 
nitely. The upper graph of Figure 3.3 confirms this result. 
The variance of the temperature perturbation can clearly be 
seen to exhibit a ramp type growth over the period of the 
experiment. Also, there is no statistically significant 
difference between the average temperature with random 
forcing and the average temperature response to the annual 
cycle alone, T-To. If the experiment was repeated an 
infinite number of times, the values of T-To would be zero. 
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Figure 3.2 Typical year with stochastic forcing. a) 
The same annual heat flux cycle shown in Figure 3.1, 
^ random daily heat flux component Q added, 
over the model year. The mixed layer depth 

h=100 m. 
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Figure 3.3 Statistical temperature variables. a) The 
mean temperature variance ot 100 stochastically forced 
cases. b) Plot of T-T<, . 
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This result only holds for the case where the mixed layer 
depth is fixed, which is a physically unrealistic situation. 

From the above results, it can be said that the 
modified Garwood model in slab form correctly duplicates the 
linear time rate of change of SST developed by Frankignoul 
and Hasselmann (1977). Because the Garwood model is not 
linear and the forcing is completely specified, it is not 
possible to compute an adequate negative feedback term, 
without first determining the climatology, to confirm 
further studies made by Frankignoul and Hasselmann. In the 
next section, we explore the role of seasonality in the 
ocean response to stochastic forcing using the one- 
dimensional model with entrainment into the mixed layer. 

3 . Stochastic Forcing with Mixed Layer Physics 

For this experiment, the governing equations are 
(2.1), (2.2), and the rate of entrainment is given by equa- 
tion (2.4). This study presents a significant departure 
from the slab models used by Hasselmann (1976) and 
Frankignoul and Hasselmann (1977), because the mixed layer 
temperature is no longer independent of the mixed layer 
depth. Given a constant wind stress, the mixed layer depth 
will shallow during periods of strong downward surface heat 
flux, which also causes an increase in mixed layer tempera- 
ture. Thus on an annual basis, the effects of the heat flux 
and shallow mixed layer depth will reinforce each other and 
generate larger perturbations in temperature. The atmos- 
pheric forcing and model initializations will remain as 
stated for the previous study. 

Figure 3.4 displays the model response to the annual 
period atmospheric surface heat flux without a stochastic 
forcing component. As stated earlier the atmospheric 
forcing amplitude is designed to allow the SST to increase 
by approximately 8 C and the mixed layer depth to shallow to 
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Figure 3.4 Annual cycle without stochastic forcing, 
a) The daily SST response to the annual cycle. b) the 
daily mixed layer depth response to the annual cycle. 



30 



approximately 50 m. These conditions are met in this 
version of the model. 

The SST reaches a maximum at approximately day 170 
which indicates a 80 day lag from the point of maximum down- 
ward heat flux. As the second model year of each case is 
used for calculations, the mixed layer depth (Figure 3.4) at 
vernal equinox is slightly shallower than the 100 m used to 
initialize the model. The depth continues to shallow and 
reaches minimum values at about day 100. For the remainder 
of the model year, the mixed layer depth gradually deepens 
until approximately day 350. From day 310 to day 350, the 
mixed layer depth has achieved the greatest depth that 
turbulent mixing will allow, and it is entraining the top of 
the permanent thermocline. During the last 10 days, the 
mixed layer begins to shallow toward the initial value. 

Upon comparing Figures 3.4 and 3.1, it is apparent 
that the maximum SST is greater when the mixed layer depth 
is allowed to vary. As the atmosphere adds heat to the 
ocean, the mixed layer becomes more stable, and the poten- 
tial energy of the system increases. Increasing the poten- 
tial energy of the mixed layer is at the expense of the 
turbulent kinetic energy, which causes the mixed layer depth 
to retreat according to equation (2.5). The initial shal- 
lowing of the mixed layer corresponds to the period of 
reduced surface cooling in late winter and early spring when 
entrainment ceases because of the reduced buoyancy produc- 
tion of turbulent kinetic energy and a relative increase in 
the dissipation of turbulent kinetic energy. Equation (2.1) 
can be reduced, for this case with forced heat fluxes, to a 
ratio of the atmospheric heat flux and the mixed layer 
depth. As the time rate of change in temperature is 
inversely proportional to the mixed layer depth, the shal- 
lowing of the mixed layer coupled with the downward heat 
flux produces a rapid increase in temperature with time. 
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Fieure 3 5 Typical year with, stochastic 
mile! k) The SST over the model year, 
layer depth over the model year. 
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Figure 3.5 is for a typical model year under the 
influence of both seasonal and stochastic forcing. The 
daily random heat flux component Q' is identical to the 
forcing in Figure 3.2. In all of the graphs the annual 
cycle is evident. The first significant result can be seen 
in the graph of mixed layer depth for a typical year in 
Figure 3.5. It is evident that the daily variations in 
mixed layer depth can be as great as 50 m during portions of 
the year. Therefore, a linear model of stochastic forcing 
with fixed layer depth is unrealistic. The temperature 
perturbation is dependent on the interaction between the 
annual plus stochastically forced heat fluxes and the 
perturbations in the mixed layer depth. 

Figure 3.6 shows the temperature statistics calcu- 
lated from the 100 cases in the experiment. There is a 
phase shift associated with the time of maximum SST, with 
the maximum SST occurring about two weeks later for the 
stochastically forced cases. The difference between this 
result and that obtained for the slab model, can be attrib- 
uted to the combined effects of a changing mixed layer depth 
and the evolving SST. This curve does retain a strong 
annual cycle. 

The temperature variance curve of Figure 3.6 retains 
the increased amplitude throughout the year observed in 
Figure 3.3 for the slab model, but the influence of the 
annual cycle is clearly evident and produces a 1 C maximum 
temperature variance. The growth visible in the variance 
graph is a result of the lack of any negative feedback 
applied to the model. 

To further delineate the temperature differences, 
the graph of T-To reveals that the stochastically forced 
average temperature is approximately 0.75 C higher than the 
temperature derived from the annual cycle alone. This 
result is significant. Both curves have a strong annual 
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Figure 3.6 Temperature statistics-full model. a) The 
variance of the temperature perturbation, b) The dif- 
ference between the mean temperature resulting from 
the stochastically forced cases and the ocean response 
to the annual cycle To . 
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cycle, which indicates that the seasonal variations cannot 
be ignored when studying SST perturbations. 

Figure 3.7 shows the mixed layer depth statistics 
analogous to those computed for temperature. It is apparent 
that the average mixed layer depth with stochastic forcing 
is smaller than the mixed layer depth produced by the annual 
cycle alone. A considerable phase shift in the minimum 
mixed layer depth occurs. The stochastically forced depth 
minimum occurs approximately one month later than the 
minimum in h<> . It is this phase difference which produces a 
corresponding phase shift in the temperature statistic, 
T-To. As a result of this phase shift a second harmonic 
becomes evident in the h-h<» curve. The final effect of note 
pertains to the mixed layer depth variance. The mixed layer 
exhibits greater fluctuations during the shallowing regime 
than the deepening regime. This phenomenon can be explained 
through an examination of the turbulent kinetic energy 
budget. In the spring the near-zero heat flux Q©, allows 
the mixed layer depth to be controlled by Q' , the random 
heat flux component. At this time, h' will show a large 
variation in proportion to Q' . In addition, there is a 
small temperature jump (AT) across the seasonal thermocline 
and, as can be seen in equations (2.2) and (2.4), the time 
rate of change of mixed layer depth is inversely propor- 
tional to AT. In the fall, the temperature jump across the 
seasonal thermocline is very large. This greatly restricts 
mixed layer depth fluctuations. ‘ Thus the shallowing regime 
has a greater mixed layer depth variance than does the deep- 
ening regime. The curve of h-h® is always negative, except 
for a random phenomenon occurring around day 20, which indi- 
cates that the stochastically forced mixed layer is on 
average 20 m shallower than the mixed layer produced solely 
by the annual cycle. 
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Figure 3.7 .Mixed. layer depth statistics- full model, 
a; The variance in the mixed layer depth over the 
experiment. b) The difference between the average 
mixed layer depth when randomly forced and the mixed 
layer depth which is the product of the annual cycle 
alone, h-ho. 
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B. INTERPRETATIONS 



1. " slab Model " Experiment 



With a constant depth mixed layer the results 
obtained confirm those of previous researchers. The ocean 
does respond to white noise atmospheric forcing with a 
"red-shifted" spectrum. The annual cycle affects only the 
average temperature, not the variance in temperature. 
However, the temperature variance does display the indefi- 
nite growth predicted by Hasselmann, since there is no 
negative feedback. 

The assumptions used in this experiment produce a 
linear equation for the change in T’. There are no 
temperature-mixed layer depth interactions. This simple 
model eliminates the possible contribution of the applied 
annual period heat flux forcing. The results would be the 
same if the random surface heat flux component were applied 
to a constant surface heat flux. Therefore, the initial 
model assumptions of Hasselmann and Frankignoul and 
Hasselmann oversimplify the problem. To produce more real- 
istic results requires the added complexity associated with 
the nonlinear dependence between mixed layer depth and SST. 

For simplicity, assume that the time rate of change 
of SST is given as follows: 




(3.1) 



By means of a Reynolds decomposition it can be shown that it 
is possible to linearize equation (3.1) only if I h’ I << h. 
Expanding equation (3.1) to include the average and pertur- 
bation terms results in the following expression 
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(3.2) 




^ o' 

(1 ♦ ‘"/h) • 



A binomial expansion gives 



/ ./ . .2 . 

0 ♦ -S-) : 1 - ♦ (^) - •• • 1 ■ r • (3.3) 

h h K 



Therefore the approximate form of equation (3.1) becomes 



)tT . ■ Q»Q’ m , il.') (3.4) 

it * At - ;; h 



If Ih'l << h, equation (3.4) reduces to 




(3.5) 



Averaging each term produces 




(3.6) 



Subtracting (3.6) from (3.5) results in the following linear 
equation for the change in T' 



This equation is the functional form used by Hasselmann 
(1976). It is evident that the temperature perturbation 
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(3.7) 



- _Q_ 
■ h 



will be the integral of the right hand side. Thus the ocean 
acts as an integrator of the random forcing. It is equally 
evident that the annual heat flux cycle plays no part in the 
determination of the temperature perturbation. 

2 . One -Dimensional Model Experiment 

Suppose Ih'l is not much less then h. Expanding 
equation (3.4) gives 



= Q . Qk .. - Q b 



Averaging each term produces 




Subtracting (3.9) from (3.8) gives 



-Q1 . JUl . i2V . iiJl' 



(3.8) 



(3.9) 



(3.10) 



Here the time rate of change of T' depends on the correla- 
tion between the random forcing component and the perturba- 
tion mixed layer depth and on the relationship between the 
average heat flux and the mixed layer depth perturbation. 
If, as in the second experiment, Ih'l is not much smaller 
than h, a linear relationship for T' is not possible. 
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The results of these two experiments clearly demon- 
strate the complex relationships existing in a single phase 
of ocean response. The linear approach to stochastic 
forcing does bring out some basic responses, but ignores 
some very significant interactions. The ocean does respond 
as an integrator of the random forcing and the variance in 
SST does grow indefinitely when no feedback is applied, but 
the strong influence of the annual cycle is missed when a 
slab model is used. In addition, the dynamics of the slab 
model are such that the average of the stochastically forced 
climatology is equivalent to the annual cycle climatology. 
This is not realistic. The two climatologies are not iden- 
tical. The random forcing component imparts enough extra 
energy to the ocean mixed layer to raise the temperature by 
0.75 C during the summer, and it delays the temperature 
maximum and depth minimum until later in the year. These 
effects are important. 

C. STOCHASTIC WIND STRESS EXPERIMENT 

The next step in the investigation of the stochastic 
forcing theory is to examine the model response to a wind 
stress composed of a mean component and a random perturba- 
tion. Previous investigations into the importance of the 
wind stress in the development of SST anomalies indicate 
this is an important forcing mechanism (Camp and Elsberry, 
1978; Frankignoul and Reynolds, 1983). 

1 . General Method 

A modified version of the Garwood (1977) model 
described in the previous sections is used for this investi- 
gation. The model is modified such that there is no longer 
a stochastic heat flux component Q' , but the annual period 
heat flux cycle remains. The kinematic wind stress is 



40 



decomposed into a mean component U** and a stochastic compo- 
nent (U*^)'. The mean kinematic wind stress is 1.0 cm*/s^ 
and the random component is allowed a daily variation of 
plus or minus 0.9 cm^/s^. In all other respects the proce- 
dure is the same as that used in the stochastic heat flux 
investigations . 

2. Stochastic Wind Stress with Seasonal Thermo dine 

This experiment isolates the model response to 
stochastically forced wind stress and a simple sinusoidal 
annual surface heat flux while retaining a variable mixed 
layer depth. As in previous investigations, the first case 
possesses no stochastic forcing to determine the base or 
zero state. The model is then run 100 years and the output 
averaged to approximate climatology. The temperature vari- 
ance (Figure 3.8) does not appear to display an indefinite 
growth in spite of the fact that there is no feedback 
applied to the model. Such growth, although still present, 
is insignificant at the scale of the figure. Additionally, 
the maximum temperature variance occurs about day 112, 
approximately 40 days prior to the maximum average SST. The 
magnitude of the temperature variance is much smaller than 
for the stochastic heat flux experiment. The annual cycle 
is evident in both curves of Figure 3.8. The mean tempera- 
ture response to stochastic wind forcing is smaller than the 
temperature response to the annual cycle alone. It can be 
seen that the temperature of the stochastically forced cases 
is about 0.25 C smaller than To, the temperature response to 
the annual cycle alone. 

Figure 3.9 shows the mixed layer depth statistics 
for the random wind stress experiment. The magnitude of the 
variance of the mixed layer depth in the shallowing regime 
is greater than it was for the heat flux experiment (Figure 
3.7). During the deepening regime, the variance in h'* 
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3 . 9 Stochastic winds ~mixed layer depth statis- 
tics, a) The variance in the mixed layer depth over a 
100 year climatology b) The difference between the 
mixed layer depth when randomly forced and the mixed 
layer depth response to the annual cycle? mixea 
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approaches zero. This phenomenon can be explained through 
an examination of the turbulent kinetic energy budget (see 
discussion for Figure 3.7). Examination of the lower graph 
of Figure 3.9 reveals that the stochastically forced mixed 
layer depth is about 15 m smaller than h during the shal- 
lowing regime. In comparison with Figure 3.7, it is evident 
that the mean mixed layer depth with stochastic wind stress 
forcing is about 5 m smaller than with stochastic heat flux 
forcing. During the deepening regime, however, the mean 
mixed layer with stochastic wind forcing is actually deeper 
than h the model response to the annual cycle. Therefore, 
the stochastic wind stress component causes the mixed layer 
to deepen more rapidly through additional generation of 
turbulent kinetic energy that is not completely compensated 
by dissipation. 

In conclusion, it can be stated that the stochastic 
forcing concept when applied to the wind stress U** does 
affect the SST response of the ocean model by as much as 
0.4 C. The stochastic wind stress forcing reduces SST on 
the whole because of enhanced mixing. The random wind 
stress forcing becomes an important source of turbulent 
kinetic energy during the deepening regime in the late 
summer and early winter. The mixed layer deepens more 
rapidly and achieves a greater depth than without the 
stochastic wind stress component. 



44 



IV. SUMMARY AND RECOMMENDATIONS 



The Hasselmann stochastic forcing theory was verified 
for a slab model, and then it was modified to include a 
mixed layer depth model. This enabled the exploration of 
the interaction between stochastic and annual period forcing 
of the ocean. The simple linear prediction equation for T' 
that was developed by Frankignoul and Hasselmann was no 
longer applicable because of nonlinear interactions between 
h’ and T'. The mixed layer depth perturbations are signifi- 
cant in comparison with the average mixed layer depth. The 
temperature perturbation is a function not only of the 
stochastic heat flux component, but also of the interaction 
between the annual heat flux cycle Q* , the stochastic 
forcing component Q' , and the mixed layer depth perturba- 
tion. Of considerable significance is the fact that with 
the stochastic forcing the climatological SST response is 
0.75 C higher than the corresponding ocean response to the 
annual cycle alone. This is consistent with the climatolo- 
gical mixed layer depth being smaller when random forcing is 
applied. Furthermore, the minimum mixed layer depth occurs 
almost one month later in the year than the ocean response 
to the annual cycle alone. Therefore, the stochastic 
forcing component of the atmospheric heat flux greatly 
influences the climatology of both the SST and mixed layer 
depth. 

In confirmation of Elsberry and Camp (1978), we show 
that the wind stress with a stochastic component produces a 
cooler (by 0.25 C) ocean surface than would be observed 
without the random component. This is because the mixed 
layer is deeper on the average due to the generation of 
additional turbulent kinetic energy. Interestingly, however. 
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the climatological mixed layer depth is reduced in value 
during much of the cooling season. 

Recommended future effort will be to apply the 
stochastic forcing simultaneously to the wind stress and the 
heat flux. This represents the most realistic situation, 
but care must be exercised to insure the proper correlation 
between the random forcing of each variable, because wind 
stress and heat flux are not independent. This is the first 
step to analyzing the effect of stochastic forcing on the 
individual turbulence components. 

The SST response of the ocean on climatological time 
scales is dependent on both annual cycles and stochastic 
forcing components. The interactions are complex and 
nonlinear. Only by examining the parts is it possible to 
understand the whole, but care must be exercised to insure 
that the simplifications used do not misrepresent the actual 
interactions . 



46 



REFERENCES 



Bjerknes, J., 1966: "A Possible Response of the Atmospheric 
Hadley Circulation to Equatorial Anomalies of Ocean 
Temperature." Tellus , 28 , 820-828. 

Camp, N. T., and Elsberry, R. L. , 1978: "Oceanic Thermal 
Response to Strong Atmospheric Forcing II. The Role of 
One-Dimensional Processes. J. Phys . Oceangr . , 8, 215-224. 

Denman. K. L. , 1973: "A Time -Dependent Model of the Upper 
Ocean.’' J. Phys . Oceangr . , 3, 173-184. 

Elsberry, R. L. , and Camp, N. T., 1978: "Oceanic Thermal 
Response to Strong. Atmospheric Forcing I. Characteristics 
of Forcing Events. j. Phys . Oceangr . , 8, 206-214. 

Frankignoul, C. , and Hasselmann, K. , 1977: "Stochastic 
Climate Models, Part II Application to Sea Surface 
Temperature Anomalies and Thermocline Variability." Tellus, 
^ 14), 289-305. 

Frankignoul, C., and Reynolds, R. W. , 1983: "Testing a 
Dynamical Model for Mid-Latitude Sea Surface Temperature 
Anomalies.' J. Phys . Oceangr . , 13 , 1131-1145. 

Garwood, R. W., 1977: "An Ocean Mixed Layer Model Capable 
of Simulating Cyclic States. J. Phys. Oceangr., 7, 
455-468. - — ^ ^ 

Garwood, R. W. , 1979: "Air-Sea Interaction and the Dynamics 
of the Surface Mixed Layer.' Rev. Geophys . Space Phys., 17 
(7), 1507-1524. ^ — 

Garwood, R. W. , and Yun, J. “Y., 1979: "Bulk Closure for 
the Oceanic Planetary Boundary Layer -- A Realistic and 
Numerically Efficient Model." Proc. of 2nd Symposium on 
Turbulent Shear Flows, Imperial College~Eondon, July I?79T 
12.6-12. 11“ 

Garwood, R. W., and Adamec, D.. 1981; Naval Postgraduate 
School Technical Report NPS68-82-006 , Model Simulation of 
Seventeen Years of Mixed Layer Evolution at Ocean Station 
P' apa ~ 1^37-^ — 

Hasselmann. K. , 1976: "Stochastic Climate Models-Part I. 
Theory." tellus , 28 , 473-485. 

Kraus, E. B., and Turner, J. S., 1967: "A One-Dimensional 
Model of the Seasonal Thermocline: II. The General Theory 
and its Consequences." Tellus , 19 , 98-106. 

Lemke , P., 1977: "Stochastic Climate Models, Part 3. 
Application to Zonally Averaged Energy Models. Tellus , 29 
(5j, 385“392. 

Lilly, D. K. , 1972; "Numerical simulation Studies of Two- 
Dimensional Turbulence." Geophys . Fluid Dyn . , 3 , 289-319. 

Namias , J., 1959: "Recent Seasonal Interactions over the 
North Pacific Waters and the Overlying Atmospheric 
Circulations.’ J. Geophys . Res . , 64 , 631-646. 



47 



Namias , J., 1969: "Seasonal Interactions Between the North 
Pacific Ocean and the Atmosphere During the 1960 's.' Mon. 
Wea . Rev. . 97, 173-192. 

Rhines. P. B., 1975: "Waves and Turbulence on a Beta 
Plane.’' J. Flui^ Mech . , 69, 417-443. 



Rowntree, P. R. , 1976: "Response of the 
Tropical Atlantic Ocean Temperature Anomaly. 
Roy . Met . Soc . , 102 , 607-625. 



. Atmosphere 
" Quart . 



to 

J. 



Shukla, J., 1975: Effect of Arabian Sea-Surface 

Temperature Anomaly on Indian Summer Monsoon: 
Numerical Experiment with the GFDL Model." J. Atmos. Sci., 
503-511. - 

Wyrtki, K., 1965: "The Average Heat Balance of the North 

Pacific Ocean and its Relation to Ocean Circulation." J. 
Geophys . Res . , 70 (18), 4547-4559. ~ 



48 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, VA 22314 

2. Library, Code 0142 ' 2 

Naval Postgraduate School 

Monterey, CA 93943 

3. Professor Robert J. Renard, Code 63Rd 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943 

4. Professor Christopher N. K. Mooers , Code 68Mr 1 
Department of Oceanography 

Naval Postgraduate School 
Monterey, CA 93943 

5. Professor Russell L. Elsberry, Code 63Es 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943 

6. Mr. Patrick C. Gallacher, Code 63 1 

Department of Meteorology 

Naval Postgraduate School 
Monterey, CA 93943 

7. Professor Roland W. Garwood, Code 68Gd 2 

Department of Oceanography 

Naval Postgraduate School 
Monterey, CA 93943 

8. LT. David C. Copley 2 

U.S. Navy Oceanography Command Center 

Box 12, COMNAVMARIANAS 
FPO San Francisco, CA 96630 

9. Directory Naval Oceanography Division 1 

Naval Observatory 

34th and Massachusetts Avenue NW 
Washington, D.C. 20390 

10. Commander 1 

Naval Oceanography Command 

NSTL Station 

Bay St. Louis, MS 39522 

11. Commanding Officer 1 

Naval Oceanographic Office 

NSTL Station 

Bay St. Louis, MS 39522 

12. Commanding Officer 1 

Fleet Numerical Oceanography Center 

Monterey, CA 93943 



49 



13. 



1 



14. 

15. 

16. 

17. 

18. 

19. 

20 . 
21 . 
22 . 

23. 

24. 



Commanding Officer 

Naval Ocean Research and Development Activity 

NSTL Station 

Bay St. Louis, MS 39522 



Commanding Officer 

Naval Environmental Prediction Research 
Facility 

Monterey, CA 93940 

Chairman, Oceanography Department 
U.S. Naval Academy 
Annapolis, MD 214D2 

Chief of Naval Research 
800 N. Quincy Street 
Arlington, VA 22217 

Office of Naval Research (Code 420) 

Naval Ocean Research and Development Activity 
800 N. Quincy Street 
Arlington, VA 22217 

Scientific Liaison Office 
Office of Naval Research 
Scripps Institution of Oceanography 
La Jolla, CA 92037 



Library 

Scripps Institution of Oceanography 
P.O. Box 2367 
La Jolla, CA 92037 

Professor Robert L. Haney, Code 63Hy 
Department of Meteorology 
Naval Postgraduate School 
Monterey, CA 93943 

Library 

Department of Oceanography 
University of Washington 
Seattle, WA 98105 



Library 

School of Oceanography 
Oregon State University 
Corvallis, OR 97331 



Commander 

Oceanographic Systems Pacific 
Box 1390 

Pearl Harbor, HI 96860 

Professor Calvin R. Dunlap, Code 68Dn 
Department of Oceanography 
Naval Postgraduate School 
Monterey, CA 93943 



1 



1 



1 



1 



1 



1 



1 



1 



1 



1 



1 



50 




3 6 



3 7 f 



i 



I 



2iuib'3 



Thesis 

c.l 



Thesis 

CT543 

C.l 



Copley 

The effect of stoc- 
hastic sTjrface heat 
fluxes on the climato- 
logy of the seasonal 
thermocline. 



2.LU153 



Copley 

The effect of stoc- 
hastic surface heat 
fluxes on the climato- 
logy of the seasonal 
thermocline. 



thesC7543 

The effect of stochastic surface heat fl 




3 2768 002 09142 3 

DUDLEY KNOX LIBRARY 








